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Introduction 



The Gross-Neveu model is the simplest version of models with four-fermion interactions 
first introduced in [Q] as examples of dynamical symmetry breaking. These models have 
been widely studied in the literature in many different areas of theoretical physics: for 
example as low energy effective theories of QCD or, in the Euclidean formalism, as 
realistic models describing phase transitions in statistical mechanics (see for example [^], 
where the stability conditions for the effective action in the large N limit are seen to imply 
the BCS gap equation of classical superconductivity). In the cosmological framework the 
massive composite field, a, which arises after dynamical symmetry breaking, has been 
proposed as an inflaton in inflationary universe models 0. With respect to each of 
the mentioned problems it is of interest to consider manifolds other than i?". In facts, 
the compactification of one direction from R to allows to describe n — 1 dimensional 
systems at finite temperature. The introduction of scale parameters in the other directions 
allows, in low dimensional physics, the study of the behaviour of dynamical systems under 
deformation of their microscopic structure. Finally, the most obvious application of the 
study of this model in the non-zero curvature case is the cosmological framework. 

There exists a very reach literature on four fermion models. An introduction is con- 
tained in p. A recent and updated review may be found in |^. The issue of critical 
exponents and j3 function in flat space is addressed for example in [|, ^ |TU|. In |TT|-P5 



the effects of an external electromagnetic field in fiat space-time are considered. The proof 
of 1/iV renormalizability in 3 dimensions was first addressed in [|T6|, [1^ . The thermo- 
dynamical behaviour in flat space-time is described in |T9|-|2^, where the temperature 



induced second order phase transition is found. The influence of a classical gravitational 



field on the dynamical symmetry breaking has been analyzed in |^-|^, where the exis- 
tence of curvature induced phase transitions is evidentiated. The relevance of an external 
electromagnetic field to descriptions of the early universe is considered in a series of pa- 
pers where the electromagnetic effects are studied in combination with curvature effects 



2mm 



In this paper we study the large N limit of the GN model, on 3-d manifolds of constant, 
non-zero curvature, in the Euclidean formalism. We focus our attention on the second 
order, temperature induced, phase transition and the role of curvature in the process of 
symmetry restoration. This may be formalized by considering the model on manifolds 
of the form S x 5*^, where S is a 2-dimensional surface with non-zero curvature and the 
inverse radius of the circle plays the role of the temperature (the compactification of one 
space-time dimension to the circle is also referred as a finite size effect). To our knowledge. 



the combined effect of curvature and finite temperature have only been considered in p2 
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(in the weak curvature approximation) and [0. In |Q the GN model is considered on 
manifolds of the form S x 5*^, but the thermodynamical phase transition is not discussed 
and the coupling constant is fixed to the fiat space critical value. In non-perturbative 
approaches it is a known result, although yet debated, that positive curvature and fi- 
nite temperature enhance the process of symmetry restoration, while negative curvature 
favours symmetry breaking. For positive curvature, a finite temperature does not mod- 
ify the qualitative features of the model but simply changes the phase transition point, 
whereas, for negative curvature, there is no phase transition at T = (the symmetry is 
broken for any value of the coupling and of the curvature), but as soon as the temperature 
is switched on a symmetry restoration becomes possible, indeed it is realized, for some 
value of T. 

We consider the GN model in de Sitter and anti-de Sitter backgrounds, respectively 
the manifolds x 5*^ and x S^. We apply the large N approximation, and use 
the zeta function regularization scheme. The GN model exhibits on a two-phase 
structure, the phase transition occurring for a non-trivial value of the coupling constant, 
the ultraviolet (UV) fixed point. When considering the model in curved backgrounds, 
at finite temperature, the Lagrangian turns out to be dependent on three parameters: 
the coupling constant, q, the parameter (3 (the inverse temperature), and the curvature 
parameter, r. Phase transitions can in principle occur with respect to any of them. With 
respect to the thermodynamical characteristics, we recover the qualitative behaviour of 
[0], and we find a simple analytic expression for the critical surfaces f{r,(3,q) = 0. We 
follow the conventions of and |^3[ where the zeta function regularization is adopted 



to study the non-linear sigma model and the GN model, respectively, at the critical value 
of the coupling constant. We find some interesting common features between the GN 
model and the non-linear a model. For example, the mass gap on x R at the critical 
coupling, for the GN model is found to be half the mass gap of the non-linear a model 
on the same manifold. 

In section 1 we briefly review the properties of the GN model and the large N limit. 
Then we recall the flat space analysis in the zeta-function regularization. The curvature 
and temperature effects are considered in section 2. The two appendices are devoted to 
calculations related to section 2. 



1 The Gross— Neveu model in three dimensions 

In this section, upon giving some basic definitions for the study of the model on a generic 
Riemannian manifold in the expansion, we review the zeta function regularization 
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scheme as discussed in |Q. Then, to fix the notation and to illustrate the method in a 
simple case, we analyze the model in the Euclidean fiat space R^. 

The Gross-Neveu model on a Riemannian manifold {/iA,g) is described in terms of a 
0{n) symmetric action for a set of massless Dirac fermions. The Euclidean partition 
function in 3-dimensions in the presence of a background metric g^v{x) is given by 



Z[g]= V[^]V[ilj] exp / d'x ^ 



M 



'1.11 



where i = 1, 2, ■ ■ ■ , A^, y is the Dirac operator on M.^ and q is the coupling constant[]. 
The Dirac matrices are given in terms of the Pauli matrices cTa by the expression 



= V^^a<ya , with /i, a = 1, 2, 3 
where V^^a denote the dreibein defined by the equation 

The covariant derivative acting on a spinor field is defined as 



where r„ is the spin connection 



;i.2) 



;i.3) 



(1.4) 



;i.5) 



In even dimensions the model has a discrete chiral symmetry which prevents the addition 
of a fermion mass term, while in odd dimensions a mass term breaks space parity. 

As it is usually done, the partition function ( |1 . 1| ) is rewritten by introducing an aux- 
iliary scalar field a, such that 



Z[g] = J V[^] V[^] V[a] exp | d^x ^ 



2q 



-a 



(1.6) 



The a field has no dynamical effect, its introduction amounting to multiply the partition 
function by an overall constant. We note here that in the large A^ approximation the 
generating functional (|1.6|) describes the Nambu-Jona-Lasinio model as well, when the 
space-time dimension is greater than 2. This model differs by the GN model in the 
presence of an extra term in the action. It reads 



Snjl = Sgn - d'^x ^/g {ipa^i'iY 



M 



(1.7) 



t According to our notation the Dirac matrices obey the following algebra: 7^ = jji, {ituli^} — 2(7^,^, 
and Tr (7^) — 0. Thus, the Dirac operator is antihermitian = — y. 
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Following the same method as above, it can be seen that the new interaction gives rise to 
another auxiliary field, called vr in the literature. In the large limit, for d > 2, n can be 
reabsorbed into the definition of the a field, giving rise to the generating functional ( |1.6| ) 
(see for example [0 for details). 

The generating functional Z has to be regularized in the ultraviolet. This is usually 
done by introducing a cut-off. A, in the momentum space. By redefining the dimen- 
sional coupling constant l/g(A) as A/g'(A), the regularized partition function is formally 
rewritten as 



Z[g,A] = / ^^aM ^^aM exp / d'x ^ 



A 

2^^ 



-a 



where I^a['?/'] = H dip{k) and similarly for the other fields (we will write q without the 

|A:|<A 

prime from now on). 

As for the non linear a model, the existence of a non trivial UV fixed point shows that 
the large momentum behaviour is not given by perturbation theory above 2 dimensions, 
where the theory is asymptotically free (see for example |^). We use here the 1/A^ 
expansion, which has the property of providing a renormalizable theory for the 3-d GN 



model |16, 17, 18 



In the large A^ limit, which means N oo keeping Nq{A) fixed, the generating 
functional can be calculated using the saddle point approximation. For this purpose we 
integrate over A^ — 1 fermion fields, rescale the remaining fields ipN, i^N to \/N — 1 ipj^, 
yjN — 1 ipN, respectively, and redefine (A^ — l)g(A) as q{A). Thus we get 



VaI^n] VaI^n] Va[(t] exp {{N - l)Tr logA(y+ a)} 

A 

2^' 



Z[g,A,q{A)] = 

X exp I (A^ - 1) J d^x ^ 



1.9) 



In the limit N ^ oo the dominating contribution to the functional integral comes from the 
extremals of the action. For an arbitrary metric g^jj{x)^ these are obtained by extremizing 
the action with respect to iPn{x) keeping a{x) and iPn{x) fixed and vice-versa. Hence, a 
set of equations {gap equations) is obtained 



Ga{x,x; a,g) +'ipN'4'N - 



A 



N 



a 










(1.10) 
(1.11) 

(1.12) 
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where Ga(x, x; a, g) = (x|(y+(T) ^|x)a is the two-points correlation function of the ■i/'-field, 
evaluated for x — > x'. At the saddle point Eq. ( |1.9| ) reads 

A 



Trlog^(y+a 



Z[g,A,q{A)]=exp^{N-l 
Then, the free energy density, 

W[g,A,qiA)] 



2q Jm 
log 2 



;i.i3) 



;i.i4) 



may be calculated. In the following analysis we will look for uniform solutions of the gap 
equations 

{a)=m , {^N) = b , {tPN) = b . (1.15) 

The quantities b and b represent the vacuum expectation value (v.e.v.) of the fermion 
fields, while m, if positive, can be regarded as the mass of the field fluctuations around 
the vacuum. 

Substituting these values into Eq. ( p,.14D we have 

A 



W[m,g,A,q{A)] 



\ogj^ det(y+ m) 



2g(A) 



m 



;i.i6) 



Phase transitions occur at the extrema of the effective potential, V{a), which is defined 
in the usual way: introducing a source J for the a field in the generating functional ( |1.6| ), 
we have 

Z[J] = J V[ij] V[ij] V[a] exp {-S{^, i>,a) + Ja} ; (1.17) 
we then introduce the Legendre transform 

r[a] = -logZ[J] + Ja 

where a = ^^^j^f^- The effective potential is defined to be 



:i.i8) 



V{a) 



;i.i9) 



When evaluating the partition function in the large N limit it is found that a = m and the 
effective potential is equal to minus the free energy density . Then, the gap equations, 
(|1.10| )-( P^.12[ ) are the extrema of the effective potential as well. 

Instead than using the cutoff regularization introduced above, we use the zeta-function 
regularization, which seems to us easier to handle in the presence of curvature. Following 
[P^ , the equal point Green's function G(x, x. A) for the operator y + m is seen to be 
regularized as 

Gs{x^x]m^ g) = m{x\{/S.ii2 + rn^)~'^\x) = m C,{s^x) , (1-20) 
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so that 



G{x, x; m,g) = m hm C(s, x) 

s — >1 

C{s,x) is the local zeta function 

as,x) = Y^ixi+m'r\^, 



XI 



'1.211 



;i.22) 



+ are the eigenvalues of the operator (A1/2 + '^^) and {ipni^)} is an orthonormal 
basis of eigenvectors (A1/2 is the squared Dirac operator). The sum over the eigenvalues 
includes degeneracy and in case of a continuous spectrum the sum is replaced by an 
integral. On homogeneous spaces such as the ones we will be considering in this paper, 
C{s,x) turns out to be independent of x. Moreover the equal point Green's function 
happens to be finite for the 3 will be clear in all the situations considered below 

(this is not the case in d ^ 3). 

The gap equations ( |1.10| )-( P^.12D become, in this regularization scheme 



m lim 



6(7% -m) = 
(7^r^ + m)6 = 

{^-C(.,m)}-66 = 



;i.23) 
;i.24) 

;i.25) 



where the regularized coupling A/g(A) has been replaced by l/q{s). The free energy 
density ( |1.16| ) can be written in terms of the zeta function in its turn, recalling that 



log det(y+ m) 



logdet(Ai/2 + m2) = -- 



We have then 



W[g,q,m] 



-l-C{s,m] 
as 



s=0 



2q 



d_ 

ds 



m 



s=0 



;i.26) 



;i.27) 



where 1/q = lim^^i l/q{s) is the renormalized coupling, defined by Eq. ( |1.25| ). 

Let us see how this scheme applies to the case of R^. Upon substituting the appropriate 
eigenvalues of the Dirac operator for this space, the gap equations (|L23|)-(|L25|) now read 



m lim>^ 



d^k 



mb 
1 



(27r)3 (F + m 



2\s 



mb = 
-bb 







;i.28) 
;i.29) 



Using the Mellin transform to analytically continue the zeta function, 

1 



C(^) 



J dtf ^hR3{t]x = x')e 



;i.3o) 
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where hpisit^x = x') = (Ant) is the equal point heat kernel of the spin-i Laplacian, 
and observing that b is zero, be m zero or not, we are left with 

mlim — r^T^r-T^ = 0- (1-31) 



\qis) (47r)2 r(s) 

Performing the limit we observe that m = for positive values of the coupling, while it 
may be non zero for negative values of the coupling. At m = the second derivative 
of the effective potential changes its sign, being positive when 1/q < 0, negative when 
1/q > (we recall that Eq. ( 1.31 ) is the first derivative of the effective potential, up to a 
minus sign). Thus, the system undergoes a second order phase transition and the critical 
value of the coupling is 

— = 0. (1.32) 

qcr 

The mass gap in the broken phase is given by 

Air 

mo = (1.33) 

q 

where the index zero will serve as future reference to distinguish flat space quantities. 
Recalling that the critical value of the coupling is independent on the curvature but 
dependent on the regularization, the result (|1.32|) will be valid, in our regularization 
scheme, for all the manifolds considered, while Eq. (|1.33|) will be taken as the definition 
of the renormalized coupling when it is negative. From Eq. ( 1.27 ) the free energy density 
is now easily calculated to be 

W = -^(^+'-) (1.34) 
2 Vevr qj ^ ' 

which is seen to be zero for 1/g > (m being zero), positive for \/q negative; in the latter 
case we have 

Tfi 

Wo = — . (1.35) 



2 The Gross-Neveu model in curved space-time 

In this section we study the large limit of the Gross-Neveu model on manifolds of 
the type S x S"^ where E is a two dimensional manifold of constant curvature. The 
component can be regarded either as a compact space-time dimension, or as a way of 
introducing the temperature through the inverse radius of the circle. The former point of 
view is adopted when studying finite size effects. Here we will adopt the latter point of 
view. 
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2.1 The manifold x 

This is a positive curvature space-time, of scalar curvature 71 = {r is the radius of 
the sphere). As we aheady mentioned in the introduction, it is expected that both the 
positive curvature and the finite temperature favour the symmetry restoration. This 
result is confirmed in our approach, where we find a second order phase transition at 
some finite value of T = /(r). When either the curvature or the temperature tend to 
zero, the two-phase structure persists recovering some known results. 

We parameterize this space by x'^ = (r, 6), where < r < 2tt, —tt/2 < x ^ ^/2, 
and < ^ < 2iT. The metric tensor is then defined as 



g^i.dx'' ® dx" = cos^ x dr ^ dr + r^dx ®dx + 0^d9 ® dO 
The Dirac matrices are given in terms of the Pauli matrices a a-, by 

7i = r cos X o"! , 72 = ra2 , 73 = /5o-3 ■ 

The spin connection (|1.5| ) results to be 



(2.1) 



(2.2) 



r/.(a;) = --o-gsinx 



(2.3) 



Due to the form of the spin connection, it was seen in ||33| that the vacuum expectation 
value of the ip fields has to be zero, for the first two gap equations to be satisfied. Hence 
we are left with the gap equation 



m lim|— — — C(s, ra) 
s^i\q{s) ^ ' 







(2.4) 



Observing that the heat kernel of the Laplacian on a product manifold is just the product 
of the heat kernels on the factor spaces, and observing that the equal points heat kernel 
Laplacians on S"^ and are respectively 



hs2{t) 



27rr 



1 ^ / 
2^ / exp — -t 



Y oo 
1^ -oo 



-A'K{n + ^] t 



(2.5) 
(2.6) 



the zeta-function is obtained by Eq. (|1.30| ) to be 

32s-l 



C{s,m) 



poo ^ 



2TT'^r'^T(s) Jo 



n=—oo 1=1 



(2.7) 
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Upon evaluating the integral and taking the limit s — 1 the gap equation (p.4|) reads 



m m? 



q 27T 



rao roc / ]^ \ 

J F J dx f cot(a;) j Ki{2mxr) 



+ i!L y (-1)" ^ ^ = exp f-n^V/2 + m2r2 ) , (2i 



where P indicates the principal value of the integral. The details of the calculation are 
given in appendix A. Let us analyze the possible solutions. We already found in that 



m = is the only solution at the critical coupling 1/qcr = 0. In facts, m is equal to zero 
for 1/q > 0, the RHS of Eq. ( p.8| ) being negative for any value of (3, r. Due to this, m 
can be non zero only for negative values of the coupling. Let us fix 1/g < 0. On deriving 
Eq. ( |2.8| ) with respect to m, it may be seen that, at m = the second derivative of the 
effective potential changes its sign when varying f3,r around some Pcryf'cr- This is where 
a second order phase transition occurs We propose here a simple expression for the 
critical surface r, g) = 0. To find the critical surface we use the fact that in a second 
order phase transition the mass gap smoothly disappears as approaching criticality. Thus, 
we perform the limit m — * for the equation 

mo ^ f°° , f 1 



/ P / dx (cot (a;) ) Ki{2mxr) 

Jo Jo V xj 



^ + — V (-1)" . ^ = exp (-n-VP + m^A , (2.9) 
47r 2nr ' + m?r^ ^\ r )' ^ ' 



which is obtained by Eq. (p.8|) by factorizing out a common m and replacing the negative 
coupling with its renormalized value (|1.33|) . Recalling that {x) ""^^ 2''~^{v-l)\x-'', we 



obtain 



mo = - F - (cot(x)--)+- 5:(-irexp 
irr Jo Jo X \ X J r 



nl=\ 




1 2 ^ (-1)" 

^ n=\e^ - 1 

The principal value of the integral may be performed, yielding 

J 1 (-1) 



^- = 2-2E^^;r^> (2-11) 



where the critical radii have been rescaled by the flat space mass gap, mo. The series is 
convergent, indeed its sum is a hypergeometrical function. We recall that this expression 
is only valid for negative values of the coupling. In Fig. 1 we plot the critical parameter 
'\-lfi'cr versus the critical curvature parameter The relation (|2.11|) is particularly 
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easy to handle when the zero curvature hmit and the zero temperature hmit are to be 



performed. For r — > cxo we obtain the known result [IS] 



For /3 — > oo we get 



21og2' 
1 



r 



crJs'^xR 



2mn ' 



(2.12) 



(2.13) 



reproducing the result of |26 



To evaluate the free energy density, defined in Eq. ( [1.271 ), "we need to calculate the 
derivative of the zeta-function (|2.7|) . This is done in appendix A in some detail. The 
result is 



>V(m) 



m 



m 



m 



2q Air'^r Jo 

3 



dx 



X 



cot(x) 



i X 

- + ^ 

X 3 



K2{2mxr) 



m 



1 



Utt 487rr2 2TTf3r^' 



°° f-lV 
E ^/exp 



n,l=l 



n 



p 



(2.14) 



Different limits can be performed easily (zero temperature and/or zero curvature, zero 
mass), reproducing some known results. In the limit of zero curvature and temperature 
we recover the flat space limit ( |1.34D . Here we just consider the zero mass limit, which 
allows us to calculate the free energy at criticality. For m — > 0, Eq. ( p.l4| ) simplifies 
considerably. We obtain 



Ws^y^s^m = 0) = -- 



1 



\Til3r 



n=l 



— ^cosech^(— n). 
n 2r 



(2.15) 



This result, already found in |Q, gives the free energy at the critical line, when (3 and r 
satisfy Eq. (|2.11|) . Moreover, it reproduces the correct result at zero curvature |l^, |3^ 

3 



W{m = 0)^2^,51 
where z is the Riemann zeta function. 



:^(3) 



(2.16) 



2.2 The manifold H'^ x 

We now consider the product manifold H"^ x S*^, where Hf. is the 2-dimensional pseudo- 
sphere. This is a space with constant negative curvature. Previous calculations suggest 
that the negative curvature favours symmetry breaking so that we have to expect a com- 
petitive effect between curvature and temperature [0. This behaviour is indeed confirmed 
in our analysis. We find a second order phase transition for some finite T = /(r), while 
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the symmetry is always broken at zero T. The latter case, which may be formalized by 
considering the model on the manifold x i?, is studied in some detail and the mass 
gap is evaluated. 

We parameterize H"^ as if^ = {z = {x,y), x E R, < y < 00} , while the 
circle 5*^ is parameterized as before by 6', < < 27r. The scalar curvature of is 
TZ = — where r is a constant positive parameter. The metric tensor on the whole 
manifold is then given by 



gf,„dx^ ^dx" = ^ {dx ® dx + dy ® dy^ + 13"^ dO ® d9 



(2.17) 



where x^ = {x, y, 9) with /i = 1, 2, 3. The Dirac matrices on x 5^ are given in terms 
of flat Dirac matrices by 



71 = -'^1 ) 72 = -0"2 , 73 = /?0-3 

y y 



while the spin connection is 



(2.18) 



(2.19) 



The v.e.v. of the ip fields being zero for this case too, we are left with the gap equation 



= mlimj-^ - C(s,m)| . 
The equal-points heat kernel of the spin-| Laplacian on is [p5 



(2.20) 



{t; z = z) 



2{7ft)y^ Jo 



dx X cothxexp 



— x^r^ 



(2.21) 



whereas the equal-points heat kernel of the scalar Laplacian on the circle is given by Eq. 
(|2.6|) . Substituting in Eq. ( |1.30|) , the zeta-function reads 



C(s,m) 



27r3/2/?r(s; 

00 

X exp 

n=—oo 



dt t"~^^Hj^ dx X coth 



3^ 



X 



n + - I + m ] t 



(2.22) 



Performing the integral in t and taking the limit s — > 1 the gap equation ( |2.20|) becomes: 

ml — / — (x cothx — 1) exp 
[it (3 Jo X 



-2xr\ 



/47r2 1 
-^(n + -)2 + m2 



m 



S-^>o.(l+e-)-l 



0. 



(2.23) 
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For details of the calculation we refer to appendix B. Let us discuss the solutions for (3 
finite first. The first term in Eq. ( p.23| ) is positive definite, hence, a non zero value of 
the mass is possible in principle in the whole space of the parameters r,(3,^, even at the 
critical coupling ^ = (in |]33| this point is missed), m = is a solution as well. In this 



limit the gap equation (|2.23|) becomes indeed 



{1 f°° dx 27rr 1 1 1 
- — - / — (x cothx — l)cosech(— — x) 7:log2 \=0 (2.24) 
An (3 J-oo x^ ' ^ (3 ' 27Tf3 ^ qj ^ ' 

where the factor multiplying m is finite. In facts by the study of the second derivative 
of the effective potential (minus the first derivative of Eq. ( |2.23| )), the value m = is 
seen to be either a min or a max, depending on the values of the parameters, f3,r, ^. 
This is the second order phase transition point. Let us try a more quantitative analysis. 
To find the critical surface we proceed as in the positive curvature case. The transition 
being of second order, the mass gap, defined by Eq. ( p.23| ) up to a factor of m, smoothly 
disappears as we approach criticality. Hence, performing the limit m — 0, we get 

2nP 1 f°° dx , . , ,27rr , , ^ 

= - — [x cothx — l)cosech(— — x) — log 2 (2.25) 

q 2 j-oo X jj 



This equation, which is the analogue of Eq. (|2.11|) for the negative curvature case, is 



numerically solved in the same way. Moreover, it can be solved analytically for specific 
values of the parameters. For 27rr = (3, the integral in Eq. ( p.25| ) simplifies yielding log 2. 
Hence, when the coupling is zero (the fiat space critical value), we have 



When the coupling is positive, 
whereas, for negative coupling 



£ > h '^■^«) 

Differently from what happens in the positive curvature case, we have a phase transition 
also for positive and zero values of the coupling. The latter result in particular, tells us 
that the flat space critical coupling, 1/q = 0, which is a second order phase transition 
point at zero curvature and temperature, maintains this property when r and f3 meet the 
condition (|2.26|). 



In figure 2 we plot the temperature l/Pcr as a function of the inverse radius l/vcr for 
positive, negative, and zero values of the inverse coupling. For 1/q non-zero, the two radii 
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are rescaled by the common factor 47r/|g|.The behaviour predicted by Eqs. ( p.26| )-( p.28| ), 
can be verified to hold. The zero curvature hmit of Eq. ( 2.25|) may be performed yielding 



^ ■log2 = -i (2.29) 



27r/? ° q 

which coincides with Eq. ( p.l2| ) when the coupling constant is negative, whereas it has 
no solution in the other cases. 

Let us now consider the zero temperature limit. Since the characteristics are quite 



different, it deserves a careful analysis. If we take the limit /5 — oo of Eq. (|2.24|) , we 
get a divergent result. This means that m = is never a solution of the gap equation 
at zero temperature. In other words, the symmetry is always broken, for any value of 
the coupling and of the curvature parameter (this is to be compared with the positive 
curvature case, where at zero temperature there is a curvature-induced phase transition). 
The gap equation for this case is to be obtained by Eq. (|2.23|) taking the hmit /5 — oo, 
or, more directly, by the zeta function for the space x R. This reads 



2^2 



V /"OO poo I XT \ 

Cis) H2xRis,m) = — - — -— / dt / dx X cothxexp rn^t . (2.30) 

47r"T(s) Jo JO \ ^ / 

and it is obtained by replacing the equal points heat kernel of the scalar Laplacian on 
with the one on R, in the equation ( p.22| ). Replacing Eq. ( p.30|) in the gap equation 



( |2.2(J ) and performing the integral in t we easily arrive at the gap equation 

{m r'^ dx , , , , .mil , , 

2^^L ^^^'^^-^)Ki{'^^^rn)----^=Q. (2.31) 

The integral on the LHS being positive for any value of r, and divergent for m = 0, this 
equation yields a non zero value of the mass whatever are the sign and the value of the 
coupling constant. 

The value of the mass gap is plotted in Fig. 3 as a function of the curvature parameter 
r, for different values of the inverse coupling. At the critical coupling it takes the simple 
expression 

m{q,r) ^ ^. (2.32) 
4r 

This equation, which has the right behaviour for r ^ oo, has an interesting property: it 
states that the mass gap at the critical coupling is half the mass gap which is found in 
34] for the conformal sigma model on the same space. We think that it is an interesting 
result because it is a manifestation of bosonization, although in the large approximation 
(another warning is the fact that the result (|2.32| ) is numerical, whereas the result in 
is analytic). 
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We now consider the free energy density for the case of (3 finite. The calculation is 
performed in appendix B. Taking the derivative of the zeta function (p.22| ), the free energy 
density (|1.27|) is seen to be 



W(m) 



+ 





cothx) ^ exp 



-2xr 



487rr2 12tx 487r/5r' 



log (l + 



n=0 



m 




^2.33) 



This expression, which looks complicated at a first sight, reproduces correctly some known 
limits, as can be easily checked, like, for example, the flat space limit (|1.34|) . In particular 
it allows to calculate easily the free energy density on the critical surface. Upon taking 
the limit m — > we arrive at 



+ 



3Cfi(3) 



°° dx 

~~2 



+ 



X 

log 2 



1 X , 

— I coth X 

X 3 



27rr / 27rrx\ 

1 H —X coth 



/3 



V P J 



sinh(27rrx//3) 
(2.34) 



which yields the critical value of the energy density when r , /?, and the coupling constant 
meet the conditions ( p.26|) -( p.28|) . The analogous result reported in |Q contains an error 



in the calculation. The limit of zero curvature yields the same expression as Eq. ( |2.16| ). 
The zero temperature energy is to be derived either from Eq. ( |2.33| ) in the limit (3 oo, 
or directly from the zeta function ( |2.3CI| ). By means of the relation ( |A.10 ), which is true 
in this case as well, we find 



m 



xR 



°° dx 

X 



, 1 X 

coth X 

X 3 



K2{2xrm) 



m 



+ 



m 



127r 487rr2 



f . (2.35) 
2q 



3 Conclusions 



Using the large approximation and the zeta-function regularization to study the Gross 
Neveu model in 3 dimensions, we have found that it describes a system undergoing a 
second order phase transition when considered on spaces of constant curvature, at finite 
temperature. We have found the critical surface f{r,P,q) = 0, for de Sitter and anti-de 
Sitter spaces, where r, /3, q are respectively, the curvature parameter, the inverse tempera- 
ture, and the coupling constant. In both cases (positive and negative curvature) this turns 
out to be a critical line, f{r/q,P/q) = 0, the coupling constant appearing as a common 
scale factor when finite. In the case of positive curvature this is represented by Eq. ( |2.11| ), 
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and phase transitions are allowed only in the negative region of the coupling constant. 
For the negative curvature case we have Eq. ( p.25 ) which describes phase transitions for 
an arbitrary value of the coupling. In particular at the flat space critical coupling we have 
found an explicit solution, Eq. ( p.26| ), which states that the product of temperature and 
curvature is constant on the critical line. At zero temperature the symmetry is broken 
for any value of the curvature and of the coupling. In this case it is possible to calculate 
the mass gap. We have evaluated it at the fiat space critical value of the coupling, Eq. 
(|2.32|) , and we have found that this is half the mass gap of the conformal sigma model 
when studied under the same assumptions. 

Appendix A 

This appendix contains a proof of the results ( p.8| ) and ( p.l4| ) of the subsection 2.1. 
Let us start with the definition of the zeta-function (12I1 



C(s,m) 



dt t 



;=-oo 1=1 



47r' 



n 



1x2 

2 



1) + 



(A.l) 

To exchange the sum over n with the integral in t we use the Poisson sum formula for the 
heat kernel of the spin-| Laplacian on the circle. 



^exp 



-47r^ n + 



1 + 2^(-l)"exp 



n 

"It 



so that 
C(s,m) 



f3 



+ 2£(-l)M dtf-^I^Y.leM- 

n=l 1=1 [ 



/^^ 72 I 2 n2 

—I +m p 



72/^^ I 2 n2 

I ^ + m p 



n 



(A.2) 

(A.3) 

A + B 



The first integral. A, is performed by means of the Poisson sum formula for the heat 
kernel of the spin-i Laplacian on the sphere 



^/exp 



1=1 



^ 2 



^3 f—3/2 



/oo ( 
dx X cot X exp ^ — a^^-^ j • 



(A.4) 



This is a specialization of the usual Poisson sum formula. A derivation may be found for 



example in Appendix B of pS]. This formula allows us to exchange the sum over / with 
the integral in t, upon extracting the spurious divergences which eventually show up. We 
obtain 



A 



+ 



27r2r(s) 



s-1 



cotx ) Ks_2{2mxr) 

x; 



87r3/2r(s) 



(A.5) 
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where K^{x) is the modified Bessel function. The second integral is easier to calculate. 
It yields 



^5-1/201/2-5 oo 7 I IP 



where we have used 



^2)5/2-1/4- 



(A.6) 



Jo a 



(A.7) 



The zeta-function is then rewritten as 
C(s,m) 



27r2r(s) 



/■oo /"OO / \ 

J F J dx x^~^ ^cot X j Ks_2{2mxr) 



+ 



m 



-2s+3Y^g _ 3-j ^s-l/22l/2-s 



87r3/2r(s) 



+ 



X ^ (— l)"n 

n,i=l 



n„s-l/2_ 



7r3/2r2r(s) 



l± + ^2)s/2-l/. 



Substituting in ( p.4|) and taking the limit s — 1, we get (|2.8| ). 

To obtain the expression of the free energy density ( 2.14|) , we first rewrite the zeta- 
function ( |A.1D in the form 

1 z"" 

X 

I 



C(s,m) 



+ 



+ 



27r2r(s) 

^s-l/22l/2 



P 

JO 



dx x^ ^ ( cot fx 



+ Ks-2i'2mxr) 



^3/2,2r(,) 



^ (_l)"^-l/2. 



'/2 



P + rn2)«/2-l/4 



zi7ii^.-i/2 + 



m 



-2s 



1 r(s 



8vr3/2r(s) 487r3/2 T{s) ' 



(A.9) 



This expression differs by ( [A.8|) by terms which cancel out when s ^ 1. It is obtained 
exactly in the same way as the previous one. The extra-terms are needed to regularize 
the expression when s ^ 0. When taking the derivative of this expression with respect 
to s, the only contribution at s = is given by the terms which contain the derivative of 
l/r(s), all the others being zero (they are multiplied by a common l/r(s)). To be more 
precise, 



^as)\s=o 

ds 



where we have used 



r(.)c(s) 



' d 1 



ns). 



=o = (r(s)C(s))|,=o, 



As), 



(A.IO) 
(A.ll) 
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Hence, replacing ( [A.10| ) in the definition of the free energy, ( |1.27| ), the equation ( p.l4| ) is 
obtained. 

Appendix B 

This appendix contains a proof of the equations ( p.23| ) and ( p.33| ) of subsection 2.2. 
Let us start form the zeta- function (12.221), 



C(s,m) 



27r3/2/?r(s) Jo 







oo 

X ^ exp 

n=—OD 



t 



dt f < / dx X coth X 
2 \ 

2. 



n 



(3' 



This is rewritten as 
C(s,m) -- 



X 



27r3/2/3r(s; 

oo 

exp 



/•oo f /-oo 

j dt f'^/'^l j dx{x cothx - 1) 



n + -] + m ] t 



I roo °o r / 

+ . X / t'"^ y exp - 



n + - I + m I t 



(B.2) 



To evaluate the contribution D we use the Poisson sum formula ( A. 2 ). We have 
D = i r dtt'-^'^e-"^^^ 



(47r)3/2r(^ 



(B.3) 



We then use the result ( |A.7| ) to perform the last integral (the first is just a F function). To 
evaluate the contribution C we use the result ( |A.7| ) to perform the integral in t. Summing 
up we obtain 



C(s,m) 



2r 



7r3/2/5r(s; 



POO 

/ dxix cothx — 1) 
Jo 



(B.4) 



oo 

X E 



-3/2 



/47r' 



+ 



^0 V'^(n + l)2 + m2 
1 r(s-3/2) 



s-3/2 



Replacing this expression into the gap equation ( p^.20| ) and taking the limit s — 1 we 
obtain (|2]2|). 
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To obtain the expression of the free energy density ( p.33| ), we follow the same procedure 
as in Appendix A. We first rewrite the zeta function as 

C(.,m) = ^^,7^/^ rfx(^xcothx-l--j (B.5) 

s-3/2 



n=0 



1 r(s-3/2) 3_2, , 1 ^, ,yn(^bY-^'\^ . ^ . 



(47r)3/2 r(s) 27r3/2r(s) Y V2^. 



This expression differs from ( p.4|) by terms which vanish in the limit s — > 1. As in the 
calculation of the energy for the positive curvature case, these terms have been introduced 
to perform the limit s — > 0, without introducing divergences when exchanging integrals 
and sums. 

Observing that the derivative of the zeta-function is given by ( |A.1CI| ) and considering 
the limit s — > 0, the wanted expression, Eq. (|2.33|), is obtained. 



Acknowledgements 



The author wishes to thank M. Asorey and J. Carinena at the Department of Theoretical 
Physics of the University of Zaragoza for warm hospitality and financial support while 
this work was completed. Support from INFN is also acknowledged. 



REFERENCES 

[1] D. J. Gross and A. Neveu, Phys. Rev. DIO (1974) 3235. 
[2] Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122 (l96l) 345. 

[3] T. Muta Foundations of Quantum Chromodinamics: an introduction to perturbative 
methods in gauge theories World Scientific (1987). 

[4] L. Jacobs, Phys. Rev DIO (1974) 3956 

[5] A. M. Matheson, R. D. Ball, A. C. Davis, and R. H. Brandeberger, Nucl. Phys. 
B328 (1989) 223. 



19 



[6] J. Zinn- Justin, Quantum Field Theory and Critical Phenomena (Oxford U. P., Oxford 
1989). 



[7; 



[9 

[lo: 
[11 

[12 

[13: 

[14 
[15 
[16 
[17; 

[18' 
[19 

[2o; 

[21 
[22 
[23 
[24 

[25 
[26 



T. Inagaki, T. Muta, and S. D. Odintsov, Prog. Theor. Phys. Suppl. 127 (1997) 93 

J. A. Gracey, J. Phys. A23 (1990) L467; Phys. Rev. D50 (1994) 2840; Int. J. 
Mod. Phys. A9 (1994) 567. 

J. Zinn- Justin, Nucl. Phys. B367 (1991) 105. 

N. A. Kivel, A. S. Stepanenko and A. N. Vasilev, NucL Phys. B424 (1994) 619. 
K. G. Khmenko Z.Phys. C37 (1988) 457; ibid, C 54 (1992) 323. 
K. G. Khmenko, Theor. Math. Phys. 89 (1992) 1161. 

K. G. Khmenko, V. V. Magnitskii, and A.S. Vshivtseev, Theor. Math. Phys. 101 
(1994) 1436. 

S. Kanemura, H. Sato, and H. Tochimura, NucL Phys. B517 (1998) 567. 

E. Ehzalde, S. P. Gavrilov, S. D. Odintsov, and Yu. Shil'nov, [iiep-ph /9807368 . 

G. Parisi, Nucl. Phys. BlOO (1975) 368. 

D. J. Gross, in Methods in Field theory 1975, Les Houches Lectures, R. Bahan and 
J. Zinn- Justin eds. (North-Holland, Amsterdam, 1976). 

B. Rosenstein, B. J. Warr and S. H. Park, Phys. Rev. Let. 62 (1989) 1433. 

B. Rosenstein, B. J. Warr and S. H. Park, Phys. Rev. D39 (1989) 3088. 

G. Gat, A. Kovner, B. Rosenstein and B. J. Warr, Phys. Let. B240 (1990) 158. 

I. V. Krive and S. A. Naftulin, Phys. Rev D46 (1992) 2737. 

T. Inagaki, T. Kouno, and T. Muta, Int. J. Mod. Phys. AlO (1995) 2241. 

I. L. Buchbinder and E. N. Kirillova, Int. J. Mod. Phys. A4 (1989) 143. 

E. Elizalde, S. Leseduarte, S. Odintsov, and Yu. Shil'nov, Phys. Rev. D 53 (1996) 
1917. 

T. Inagaki, Int. J. Mod. Phys. A 11 (1996) 4561. 

K. Ishikawa, T. Inagaki, and T. Muta Mod. Phys. Lett. All (1996) 939. 



20 



[27] D. K. Kim, K. G. Klimenko, J. Phys. A 31 (1998) 5565. 
[28] D. Gitman, S. D. Odintsov, and Yu. Shil'nov, Phys. Rev. D54 (1996) 2968. 
[29] E. Elizalde, Yu. Shil'nov, and V.V. Chitov Clas. Quant. Grav. 15 (1998) 735. 
[30] T. Inagaki, S. Odintsov, and Yu. Shil'nov Int. J. Mod. Phys A 14 (1999) 481. 
[31] E. Elizalde and Yu. Shil'nov, |hep-th/9809203 



[32] E. Elizalde, S. Leseduarte, and S. Odintsov, Phys. Rev. D 49 (1994) 5551. 

S. Kanemura and H. T. Sato, Mod. Phys. Let. AlO (1995), 1777; All (1996) 785. 

[33] G. Miele and P. Vitale, Nucl.Phys. B494 (1997) 365 

[34] S. Guruswamy, S. G. Rajeev and P. Vitale, Nucl. Phys. B438 (1995) 491. 
S. Guruswamy and P. Vitale, Mod. Phys. Let. All (1996) 1047. 

[35] N. L. Balazs and A. Voros, Physics Reports 143 (1986) 109; E. D'Hoker and D. H. 
Phong, Comm. Math. Phys. 104 (1986) 537. 

[36] A. P. Prudnikov, Yu. A. Brychkov, and O. I. Marichev, Integrals and series Gordon 
and Breach Science Publishers (1992). 



21 



Figure 1: The critical temperature defined by ( |2.11| ), is plotted as a function of the critical 
curvature parameter. They are both rescaled as A = ir—, 4 = ■ 

i J b' pmo ' r' rrriQ 
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Figure 2: The critical temperature defined by (|2.25|) , is plotted as a function of the critical 
curvature parameter. The sohd line corresponds to ^ < 0. The dashed line corresponds 
to i > 0, whereas the dotted line corresponds to ^ = 0. When ^ is negative or positive, 
both the temperature and the inverse radius are rescaled to: ^ = ^^, p- = ^i. 
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Figure 3: The mass gap defined in Eq. (|2.31| ), is plotted as a function of the inverse radius. 
The sohd fine corresponds to ^ < 0. The dashed fine corresponds to ^ > 0, whereas the 
dotted fine corresponds to ^ = 0. When ^ is negative or positive, both the mass and the 



radius are rescaled to: m' 
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